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Quantum Hall bilayer phase diagram with respect to interlayer distance bears a remarkable sim¬ 
ilarity with phase diagrams of strongly correlated systems as a function of doping, with magnetic 
ordering on the one end and Fermi-liquid-like behavior on the other. Moreover, it has been sug- 
gestecP that a BCS correlated state of composite fermions with p-wave pairing may exist in the 
intermediate region. In the same region, an exact diagonalization study in the torus geometrjl^ 
pointed out the existence of state(s) with pseudospin spiraling order. Here we reconcile these two 
descriptions of the intermediate state by considering the underlying bosonic representation of the 
composite fermion paired state in the long distance limit, and by performing extensive exact diag- 
onalizations on the torus. We argue that the spiraling states belong to the manifold of degenerate 
ground state(s), and are a consequence of Bose condensation of the quasiparticles (with critical 
algebraic correlations) at non-zero momenta in the two pseudospin states. The spiraling states, 
generated in this way as spin-textures, can be identified with meron-antimeron constructions. Thus, 
merons - the fractionally charged vortex excitations of the XY magnetically ordered state - consti¬ 
tute some of the topological sectors. It follows that merons are deconfined in the intermediate state, 
and allow for a smooth transition between the magnetically ordered and Fermi-liquid-like phases, 
in which they are bound in pairs. 


I. INTRODUCTION 

Multicomponent fractional quantum Hall (FQH) sys¬ 
tems feature rich physics of strongly interacting electrons. 
In contrast to strong electronic correlations in various 
(material dependent) crystal settings, in FQH systems 
the kinetic energy is effectively removed from the prob¬ 
lem due to the Landau quantization in strong magnetic 
fields. This makes the multicomponent FQH systems a 
simpler (better controlled) setting that is amenable to the 
investigations of fundamental physics of strongly corre¬ 
lated electrons. 

One of the most prominent examples of multicompo¬ 
nent FQH systems is the quantum Hall bilayer (QHB) at 
total filling factor one. The filling factor z/ = 2tt£'^p is 
defined as a product of electron density p and the mag¬ 
netic area 2Tr£^, where £b = sJhjeB is the emergent 
length scale in the magnetic field (“magnetic length” )P. 
Alternatively, the filling factor represents the ratio of the 
number of electrons Kl to the number of magnetic flux 
quanta through the system. In the case of the v = 1 
QHB studied here, there are two quantum wells (“lay¬ 
ers”) that we label by and N = = 7V$. 

Note that throughout this work we assume the ordinary 
electron spin to be fully polarized by the external mag¬ 
netic field. 

Because of the two-fold layer degree of freedom, in 
general we need to distinguish the interaction between 
electrons in the same layer (“intra-layer”) versus elec¬ 
trons in opposite layers (“inter-layer”). The difference 
between these two interaction strengths can be tuned ex¬ 
perimentally via the parameter called the bilayer distance 
d, which makes the physics in this system much richer 
compared to an ordinary system of electrons with spin. 
QHB system has indeed attracted much attention since 


the time of the theoretical predictiorPEl and experimental 
detectiorP of the ordered state of excitons (bosons) that 
forms when d £b, despite the fact that the fundamen¬ 
tal charge carriers in the system are fermions. Perhaps 
equally remarkable and interesting is the physics at larger 
(intermediate) distances, typically 1 < d/iB ^ 2. There, 
an intermediate phase of paired composite fermions, in¬ 
troduced in Ref. [T] [see Eqs. and ([^ below], is be¬ 
lieved to exist. On the other hand, an exact diagonal¬ 
ization study of the system in the torus geometrjP also 
pointed out to the existence of pseudospin state(s) with 
spiraling order in the intermediate region. Here we rec¬ 
oncile these two descriptions of the intermediate region 
by considering the underlying bosonic representation of 
the composite fermion paired state in the long distance 
limit, and by performing extensive exact diagonalization 
studies on the torus. 

This paper is organized as follows. In Section [H] we 
start with an overview of the QHB physics from the point 
of view of model wave functions. We first discuss the 
limits of small and large bilayer distances which are well- 
understood and correspond to the exciton superfluid and 
the product of two decoupled (composite) Fermi liquid 
states, respectively. We also review the theoretical evi¬ 
dence for a paired phase in the regime of intermediate bi¬ 
layer distances, and briefly motivate and summarize our 
contributions. In Section m we provide a simple physi¬ 
cal picture for the formation of the putative intermediate 
state. Section |IV| contains results of exact diagonaliza¬ 
tion studies of the QHB on the torus. In Section jV] we 
discuss the torus results from the composite boson view¬ 
point, and give arguments for meron deconfinement in 
the intermediate region. We conclude in Section |Vl| with 
a short summary and discussion. Appendices contain 
further resullts on the intermediate phase. 
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II. BILAYER PHASES AND THEIR MODEL 
WAVE FUNCTIONS 

The QHB at the total filling factor one consists of two 
quantum Hall mono-layers separated by distance d in the 
z-direction (perpendicular to both layers). The tunneling 
between the layers can be nearly completely suppressed 
in experiment. Therefore, we can assume that = Nj^ 
and each layer is at filling one half. In the following we 
first give the basic physical picture of the ground state 
in the limit of either small {d < £b) or large distances 
d :S> £b, followed by the discussion of the intermediate 
regime 1 < < 2. 

A. Small bilayer distances: exciton superfluid 

For small distances, the intra- and inter-layer interac¬ 
tions are of similar strength, and favor each electron in 
one layer to be directly on top of a hole in the other 
layer, which drives exciton formation. In this simplified 
picture we can differentiate between up and down ex- 
citons, but due to the quantum-mechanical uncertainty 
in which layer an electron is, we must include the pseu¬ 
dospin (layer) degree of freedom. The appro pria te model 
function at total filling factor one then read^^ 

N 

'i’excizi, ■ ■ ■ ,Zn) X|-)--!■). 

i<j 

( 1 ) 

The spatial part of the wave function is written in terms 
of 2D electron coodinates Zj = Xj + iyj (regardless of the 
layer index), and assumes the circular gauge for the mag¬ 
netic field in the infinite plane. The spatial correlations 
are encoded via the Laughlin-Jastrow factor {zi — Zj), 
and the wave function vanishes at large distances because 
of the Gaussian term that originates from the single¬ 
electron orbital^. 

Note that the wave function also contains the pseu¬ 
dospin part I • • • —?►). The pseudospin of each elec¬ 
tron points to the same direction, due to the combined 
effect of the Pauli principle and the Coulomb interaction. 
The state of each particle j is given by the pseudospin 

l^.)^^(lt,)+e^'^U,)), (2) 

i.e., the pseudospin is polarized in the XY plane. This 
happens because of the capacitor effect which enforces 
the equality of total charges in the “up” and “down” 
layer, and breaks the SU(2) spin symmetry, assumed to 
exist at d = 0, down to C/(l). This correlation - condensa¬ 
tion was detected in several experiments, e.g., Spielman 
et alW, although the “zero-bias peak” anomaly was never 
as high and sharp as in the similar experiments on su¬ 
perconductors. More recent experiments using Coulomb 
drag and counterflow measurements are also completely 
in favor of exciton formation [for a recent review, see Ref. 

m- 


B. Large bilayer distances: product state of two 
composite Fermi liquids 

When the layers are far apart, it is well-known that 
each half-filled layer is in the state of the “Fermi-liquid¬ 
like” stat^. For example, the state of the upper layer is 
given by 

Nf 

■i>FL{zi,---,ZN^) = det(e5('''"^"+''’"^"^)]J(zi - Zj)^ 

i<j 

xlt-.-t), (3) 

where we dropped the single-particle Gaussian factors. 
The Fermi liquid correlations result from the factor 
det(e**^’"''") that represents the Slater determinant of 
plane waves. Here, {k} is the set of 2D momenta that 
define the 2D Fermi sphere and we have substituted 
kr = {kz + kz) /2. The total wave function of the system 
in this case is just a direct product of independent layers, 

'£!fl{z\, ..., z]^^) 0 'i'FLizt ..■,zjfj. 

Notice that both 'i>exc and dt fl contain a part which 
is the product of the differences of z’s (the Laughlin- 
Jastrow part). This part of the wave function fixes the 
filling factor of the state, which can be seen as follows. 
On the one hand, the highest exponent in any of the 
z’s is equal to the number of flux quanta through the 
system, which is the degeneracy of the lowest Landau 
level (LLL). On the other hand, this number is equal 
to the number of electrons in the case of and twice 
larger than the number of electrons in the case of as 
follows from the filling factors. Thus, formally speaking 
the Laughlin-Jastrow part is the “charge” part of the 
wave function, while the rest is the “neutral” part. In 
'kexc the neutral part is the condensate of pseudospins. 
In the neutral part contains the Slater determinant 
which introduces the complex conjugate coordinates z. 
Thus fl in general has a non-zero projection on many 
Landau levels and not just the LLL. If dtp,^ is to be used 
a LLL trial wave function, in general it must be explicitly 
projected via a complicated derivative operator Vlll to 
make it a function of z’s onljl^. Nevertheless, in the long 
distance regime with the sample length L ^ Ibi we may 
assume that z and z’s commute, and thus in this regime 
the form of the wave function formally decouples into 
charge and neutral parts. 


C. Intermediate distances: a paired phase 

Given that the QHB ground state in the limits of both 
small and large d is in a liquid phase, it seems unlikely 
that a possible intermediate state is a charge-density 
wave (that breaks the translation symmetry). 

A possible way to describe this evolution is to start 
from the Bose condensate which is smoothly modified 
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by two Fermi condensates that grow as the distance d 
increased. In such a study on the spher^ it was found 
that the two Fermi condensates pair in the manner of 
weakly paired BCS Px + ipy wave Cooper pairs. (We 
use the phrase “weakly paired” to denote Cooper pairs 
that are not exponentially localized.) Thus, at the end 
of the evolution of the superfluid state one may expect 
the following wave function, El 


N/2 N/2 

= n(4 - 4)' X n(4 - 4)' 

i<j i<i 

det(3(r^-r);)), (4) 

to describe the system. We use here and below a short¬ 
ened form of the wave functions with spin, where spa¬ 
tial coordinates carry spin. It is understood that by 
combining the spatial wave functions with spinor wave 
functions, and summing with proper sign factors over all 
ways of choosing which electrons have which spin, we 
get proper totally antisymmetric wave functions. In Eq. 
Q “det” denotes the antisymmetrized product of N/2 
pairs of Cooper pair wave functions, g{vl^ — r);), where 
m,n = 1,... ,N/2 and m ^ n. The Cooper pair wave 
function — rl^) describes Px + ipy pairing, i.e., at 

small distances it behaves aJ^^ 

5(4-4) = (^1 - ^2)'^(|ri - ral) (5) 


with /i( 0 ) 7 ^ 0 , and I = 1 which fixes the relative angular 
momentum of the Cooper pair. 

Based on the study of Ref. [TJ other numerical studies 
(summarized in Sec. |III[ ), and our results (Sec. IV), we 
will adopt the view that an intermediate state does ex¬ 
ist and can be modeled as a BCS p-wave paired state of 
composite fermions. The question of whether in the ther¬ 
modynamic limit we obtain a critical region (phase) or 
rather a critical point is more subtle because of the limi¬ 
tations of all the existing numerics to small system sizes. 
Therefore, we will use the term “intermediate phase” 
loosely, keeping in mind the possibility of it becoming a 
critical point in the limit of a very large (clean) system. 

In the remaining part of this Section, we will intro¬ 
duce in a formal way a possible alternative description of 
the intermediate state based on the picture of composite 
bosons, which we expect to be relevant in the long dis¬ 
tance limit. First, we would like to point out that the 
most natural long-distance weak pairing behavior of the 
pair function in Eq. Q with I = 1 would be 


5(4 - 4) 


1 

(zi - Z 2 ) ' 


( 6 ) 


This coincides with the existence of a BCS description 
with a gap function kx + iky, i.e., is an eigen¬ 

function of rotations in 2D with I = 1 [see Ref. E]. 

This generic behavior may be used to define a suit¬ 
able candidate for a projected p-wave paired composite 


fermion wave function valid at all distances in the LLL, 
'kpSF = T^LLLldetl^ -^ 

\ Zm Zfi 

N/2 N/2 

11 ( 4 - 4411 ( 4 - 44 }’ ( 7 ) 

i<j i<j 

where Plll denotes the necessary projection to the LLL 
(see Appendixfor the details of the projection). Here 
we confine ourselves to the long-distance regime, i.e., con¬ 
sider Eq.Q with Q. In this regime, the particles may 
be considered point-like as the relevant probing wave¬ 
lengths are much larger than £ 3 . Thus, to obtain the 
long-distance properties we use Eq. Q with ®> , instead 
of the proper LLL function [Eq. 0 ]> assuming as usual 
that the universal (topological) properties are insensitive 
to the microscopic differences between the two wave func¬ 
tions. In the following, we will explain how in this regime 
of the p-wave composite fermon pairing we can also con¬ 
sider a composite boson representation, which will en¬ 
able us to describe the topological features (degenerate 
ground states on the torus and quasiparticle content) of 
the intermediate state. 

Thus the model wave function for the state into which 
the superfluid evolves, in the long-distance limit when z 
and z commute, is 


N/2 N/2 / 1 \ 

'^psF =n (4 - 4 )^ n (4 - 4 )^^®^ ( -t _ ) ■ (^) 

i<j Kj Zn) 

Using the Cauchy identity, we can rewrite the wave func¬ 
tion in the following form 

'^PSF = n(^i - U)det ( -y—t) det ( ^ ) ■ (9) 

\Zrn- z!kJ \z'm- z/kJ 

Again, at long distances the charge and neutral part for¬ 
mally decouple, and the neutral part can be viewed as a 
wave function of some underlying composite boson quasi¬ 
particles. Furthermore, the neutral part is the product 
of two determinan ts, an d a careful conformal field the¬ 
ory (CFT) analysi^^^Esl allow for an existence of a 
branch of gapless excitations (see Appendix A) on top of 
this non-chiral state. 

These formal arguments also suggest that there is a 
possibility to interpret the intermediate phase in the 
QHB in terms of bosonic quasiparticles that correlate 
in a special way via Eq. (|^. Erom Eq. ([^ (see also Ap¬ 
pendix A) we can read off that the one-particle exciton 
(or composite boson) correlation function is 

lim (V’|(r)V’;( 0 )) ~ 4 - (Id) 

7—>-oo ' 

Therefore, instead of the long-range order that we had in 
Eq. 0, in this ground state we have quasi long-range or¬ 
der. Indeed, in the finite-size study performed in Ref. [TJ 
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it was noticed that the excitonic superfluid order param¬ 
eter, i.e. , persists in the state with no cor¬ 

related bosonic quasiparticles in the manner of Eq. Q 
(i.e., the so-called “111” stat^^^. 

It follows that these critical bosons [Eq. ([^] may exist 
in a region between a ferromagnetically ordered and a 
completely disordered phase, and thus we may expect 
a gapless mode and compressible behavior in the neutral 
sector. This is supported by exact diagonalization studies 
on the spher^i^, but also by the detection of ordered 
(correlated) behavior in early experiment^, which were 
performed at distances where we expect the physics of 
the intermediate state to be relevant. 

But what makes the intermediate region so attractive, 
as we will show in greater detail by exact diagonalization 
studies on torus [Sec. IV , is the existence of competing 
ground states at momenta that diverge in the thermo¬ 
dynamic limit, and thus define different sectors of the 
system. This signals a topological structure that co¬ 
exists with gapless behavior. A convenient way of de¬ 
scribing various sectors, as we will argue, is to take the 
composite boson view of the ground state pairing corre¬ 
lations. The n Bose condensates under different bound¬ 
ary condition j^^ l ^^ l can be identified with the competing 
ground states from exact diagonalization studies on the 
torus. 

As usual in topological systems, the torus ground 
states can be connected with quasiparticles of the topo¬ 
logical phase under consideration. We expect this also 
to be valid in our gapless system, and by Bose condens¬ 
ing in two different pseudospin states, we will be able to 
connect different ground states with meron spin-textures. 
This opens a door for a possibility that merons - frac¬ 
tionally charged vortices of the XY ordered stat^- exist 
as free (deconfined or not bound in pairs) quasiparticles 
in the intermediate state. 

The first exact diagonalization study on the torus that 
pointed out to the existence of the competing ground 
states with momenta that diverge in the thermodynamic 
limit was done by ParlP (see Appendix for a brief re¬ 
view of the relationship between Ref jU and our work). In 
Ref. [5] a spiraling structure in the competing ground 
states was identified, and it was argued that the in¬ 
termediate phase possesses a pseudospin spiral order. 
Therefore, the competing states were interpreted as new 
developing ground states that replace the XY ordered 
state (valid for smaller distances) due to a (classical) 
symmetry-breaking ordering. On the other hand, in this 
work we give numerical and analytical reasons to con¬ 
sider an intermediate state with richer topological con¬ 
tent than the XY ordered state (which has a simple in¬ 
teger quantum Hall effect characterization). In the inter¬ 
mediate phase, due to topological order, new quasiparti¬ 
cles emerge and define additional ground states on the 
torus which have a spiraling structure. In this way our 
work suggests a way to reconcile and unify the views of 
the intermediate phase based on the previous numerical 
approaches, especially Ref. [T]and Ref. [H 


III. THE FORMATION OF THE 
INTERMEDIATE STATE: A PHYSICAL 
PICTURE 


One may naively estimate that the transition to a 
disordered phase occurs around d ^ (.b, and this esti¬ 
mate is in agrement with several results in the literature 
that we summarize here, (i) The calculated treshold for 
the disappearance of the long-range correlated compos¬ 
ite bosons and the extrapolated estimate for the vanish¬ 
ing of the long-range order parameter in Ref. [TJ (ii) The 
change in the DMRG spectrum computed in Ref. [211 i-®-, 
a “change in the nature of excited state” with a “sudden 
decrease in the excitation gap”, (iii) An abrupt change 
in the velocity of the Goldstone mode in Ref. [HI (iv) Ex¬ 
act diagonalization (ED) studies (EDs) presented here as 
well as in Ref. [21 showing the level crossing for the first 
excited state around d ^ is- 

The physical picture of the evolution of the bilayer is 
the following: in the ordered XY phase we have locally 
equal probability to find an up and down exciton. As we 
increase the distance between the layers, the intra-layer 
Goulomb interaction becomes dominant and we need to 
place the up excitons, and separately down excitons, fur¬ 
ther apart. One natural resolution would be an Ising 
order (charge density wave), but another viable possibil¬ 
ity is to slowly disturb the long-range ordering and create 
spin-textures as ground state(s). In the magnetic descrip¬ 
tion of Ref. the tendency to create spin-textures may 
be recognized by the presence of frustrating next-nearest- 
neighbor Ising couplings. Thus the magnetic (bosonic) 
description of the intermediate phase may be useful along 
with the fermionic description presented in Ref. [I] 


IV. EXACT DIAGONALIZATION RESULTS ON 
THE TORUS 


In this Section we present results of exact diagonal¬ 
ization calculations performed on the torus geometry. 
We model a finite system of N electrons moving on 
a surface spanned by vectors Li and L 2 . We fix the 
filling fraction to be one, which leads to the condition 
|Li X L 2 I = 2tt£‘^N. The focus of our investigation is on 
the liquid phases, hence we consider an isotropic torus 
geometry, i.e., the length of the two vectors is the same 
(aspect ratio equal to one). Additionally, we can also 
vary the angle 9 between the two vectors Li and L 2 . 
Unless otherwise specified, we fix 0 = corresponding 
to centered rectangular symmetry. 

The interaction between electrons consists of an intra¬ 
layer Goulomb potential, given by the Fourier transform 

‘2i7T 

^intra{Q) — 

and the inter-layer Coulomb repulsion which is screened 
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FIG. 1: Energy spectrum versus momentum Mb in the odd 
sector on the torus. Panel (a) shows the spectrum for Af = 13 
electrons, while panel (b) is for = 15. In both cases the 
aspect ratio r = 1, cos0 = 0 and bilayer distance d = 1.5£b- 
Levels are plotted relative to the corresponding ground state. 


by the bilayer distance 


^inter{.Q_) 



( 12 ) 


for a recent overview of the many-body translation sym¬ 
metry classihcation in the magnetic held). 

In Figs. [2 and we present ED results on the torus 
of the bilayer system with odd numbers of particles. In 
Fig. [T] we show the energy spectrum versus conserved 
momentum Mb in the systems with Ai^ = 13 (a) and 
Af = 15 electrons (b). In this hgure, the bilayer distance 
is hxed at d = 1.5 £b, i-e., in the middle of the transition 
region. Because the total number of particles N is odd, 
we restrict to the the largest sector of the Hilbert space 
with A^.^ = {N -b l)/2 and A^^ = {N — l)/2. 

More generally, the evolution of the entire spectrum as 
the bilayer distance is varied in presented in Fig. da). 
In Fig. (b)-(d) we show several snapshots of the energy 
spectrum in Fig. [^a) plotted as a function of momen¬ 
tum. We choose three values of d = 0.5 £b, 1.5.£b, 2.5£_b 
that tentatively correspond to three different regimes: 
the 111 state, intermediate region, and the onset of de¬ 
coupled composite Fermi liquids. The unit cell in Figs. 
[^b)-(d) is approximately the square, and insets show the 
same results for the hexagonal geometry. 

What we can see in both geometries [Figs, [^b)-(d)] as 
well as Fig. is that at distance d = 1.5 £b two distinct 
states “fall” from the excited (continuum) part of the 
spectrum: for example, the state around k ^ and 
another one at A: ~ in Fig. |^c). The state with the 
momentum k ~ was first pointed out in the study of 
ParlP. Our results suggest that this state evolves hand 
in hand with another one at momentum k ^ in 

Fig. H c). The presence of these two states is robust 
and does not depend on system size, for example a larger 
N = 15 system shown in Figj^contains the same doublet 
of states. We notice that these “special” states compete 
with the ground state at k = 0 as d becomes larger. 

The two special states in fact each comprise a family 
of 4 states connected by symmetry operations (simple 
translations and inversion). Thus, in Fig. we find a set 
of 4 states that correspond to the lower energy excitation 
(with k ~ 4€^^): 

(Ar^,0), (Ar^,0), (0,iVt) and (0,iVj (13) 


The interaction potentials are made appropriately peri¬ 
odic (i.e., q assumes values 27r/Li 2 times an integer) to 
satisfy the periodic boundary conditions. 

After specifying the geometry and the interaction po¬ 
tential, the basic procedure of exact diagonalization is 
to represent the Hamiltonian of the system as a large fi¬ 
nite matrix and diagonalize it on a computer. In order 
to extract physical information about the energy spec¬ 
trum and system’s eigenstates, it is essential to label 
them with good quantum numbers. In this case, the 
complete symmetry classification was first achieved by 
Haldan^i^ the primary quantum number that emerges 
from such an analysis is the two-dimensional momentum 
that characterizes the invariance under magnetic trans¬ 
lations. This momentum is generically conserved for any 
translationally invariant interaction, and we will use it 
below to label the eigenstates of the system (see Ref. HHI 


in units of 27r/L where L = yj2TTpgN. This set of states 
is followed by the excitation at 

(iV^,iV^), (IV;,iVt)> and (iV;,IV;) (14) 

with a slightly higher energy (with momentum k ^ 6^))^). 
We note that we have identified such states in all odd 
systems we studied, e.g. N = 9, 11,13, and 15. 

Thus we may conclude that at distances where we may 
expect the intermediate phase, by examining the torus re¬ 
sults we hnd competing ground states at momenta that 
diverge in the thermodynamic limit. This may be in¬ 
dicative of the presence of Bose condensation and spin 
textures, but requires additional analysis. 

Studies of the gapless Goldstone mode, e.g. the nu¬ 
merical one on the spher^i^ and the experimental detec¬ 
tion in Ref. m suggest that even though the numerical 
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FIG. 2: (Color online) (a) Evolution of the energy spectrum in the odd sector on the torus as a function of bilayer distance 
djlB- Levels are plotted relative to the ground state at each d and quoted in units of jtin- Levels belonging to the three 
high-symmetry sectors have been designated by colours indicated in the legend, (b)-(d) Energy spectrum versus momentum 
Mb in the odd sector on the torus for d = 0.5fs, 1.5^s,2.5£s. In all cases N = 13, cos9 = 7r/2, and aspect ratio is set to 
r = 0.99 to slightly break the discrete rotation symmetry. Levels are plotted relative to the corresponding ground state for 
each value of d. Insets: Hexagonal geometry {6 = n/3 and r = 1). 



Mb 


FIG. 3: (Color online) Goldstone mode for the bilayer dis¬ 
tance d — 0.5 £b- Quasicontinuous E{k) energy spectrum is 
obtained by varying the torus angle 6 for a fixed aspect ratio 
r = 1. Data is shown for all system sizes between N = 8 and 
= 17 electrons in total. Levels are plotted relative to the 
corresponding ground state for each 6 and N. 


work, as summarized in Sec. |IV[ points to the formation 
of the disordered state around d ^ a linear gapless 
mode persists to exist in this phase. Thus although we 
may notice “a small gap structure” in Figs. after 

d ^ £b^ it is likely that it does not represent a paramag¬ 
netic phase. Because of the finite size effects and geom¬ 
etry, this “weak Goldstone” mode of, as we will argue, 
critical bosons, is distorted. Despite the distortions, in 
the hexagonal geometry at d = 2.5£b, we can see clearly, 
as we explain in the following section, consequences of 
the physics of critical bosons in the intermediate state: 
we find three (almost) degenerate ground states (families 
with members connected with symmetry operations) and 



0.02 
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Mb 

FIG. 4: (Color online) Quasicontinuous E(k) energy spec¬ 
trum obtained by varying the torus angle 6 for a fixed aspect 
ratio r = 1. Total number of electrons is = 13 and bilayer 
distance is d = 1.5^s. Levels are plotted relative to the corre¬ 
sponding ground state for each 6. Levels (A^ti 0) ^'^d (A^ti 
belong to the same excitation branch. 


discernible gapless modes that have beginnings in these 
states. 

In order to examine in more detail the question of the 
existence and description of gapless modes by EDs, en¬ 
ergy spectrum was obtained by varying the torus angle 9. 
In this way, as we sweep 9 continuously between 7r/2 and 
7 r/3, we obtain many different resolutions of momenta 
k, which produces a quasicontinuous energy spectrum. 
In Fig. we perform such a calculation in the ordered 
phase (d = 0.5£b), where we find clear evidence for a 
linear gapless mode. In Fig. |^we perform a similar cal¬ 
culation for d = 1.5£b, i.e. in the intermediate region. 
Here it is strongly suggestive that the two families of 
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states in Eqs. ( |13[ ) and ( |14[ ) indeed evolve together as 
the system is adiabatically deformed, and likely belong 
to the same gapless excitation branch. The distortions 
of a finite system do not give us a rigorous proof of the 
gapless mode, but we can expect that if such a mode ex¬ 
ists it will appear in the manner produced by the probe. 
Thus with the distortion we can just produce a state at 
the “distorted momentum”, or we can access the energy 
of the state belonging to the excitation branch in the 2D 
momentum space, at missing (undistorted) momentum 
in the finite size system. 

In the even sector {N even) of the EDs on torus, 
Fig.H we can similarly recognize the competing ground 
states with momenta that diverge in the thermodynamic 
limit. We will postpone the discussion of these results to 
Sec. VB and Sec. VC Before that, in the same section, 
it will be explained how ground states at A: 7^ 0 emerge 
from the physics of critical bosons. 
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FIG. 5: (Color online) (a) Energy spectrum versus momen¬ 
tum kis in the even sector on the torus for = 12 electrons 
in total. Aspect ratio r = 1, cos 9 = 0 and bilayer distance 
d = 1.5£b- Levels are plotted relative to the corresponding 
ground state at {N/2, N/2). Ellipses denote groups of levels 
in momentum sectors (0,0) and {N/2,0). Total degeneracy 
inside each ellipse is denoted by ID. (b) The Brillouin zone for 
= 12 corresponding to the square unit cell. 


V. A STATE OF CRITICAL BOSONS AT 
INTERMEDIATE DISTANCES 


A. The intermediate state and critical bosons 


Though from the torus results (Fig. and Fig. 
we cannot claim the ideal degeneracy of the competing 
states, their configuration is robust and characterizes the 
intermediate region. (We have verified this in systems 
of size up to 15 electrons in the odd sector, and up to 
14 electrons in the even sector.) Therefore, the systems 
analyzed are likely near a critical point, line, or perhaps 
even a regiorf^in the parameter space where there is ideal 
degeneracy, implied by the topological structure of the 
model function [Eq. ([^ or Eq. ([^]. 

In the following we explain how the “extra” competing 
ground states at non-zero momenta can be interpreted as 
pairs of merons (charged vortices of the superfluid state 
at small d), whose spiraling period is commensurate with 
the length(s) of the system. These spiraling configura¬ 
tions were considered in the work of ParlP, but here we 
would like to point out that on the torus they can be 
viewed as a consequence of the physics of critical bosons. 
In this way we will relate the physics of long spiraling 
on the torus and critical (composite) bosons to the exist¬ 
ing description in terms of p-paired states of composite 
fermions (see Introduction and Ref. [1]). 

To capture the critical correlations of bosons, in the 
following we will consider bosons with the critical BCS 
pairing function g(r) ^ l/r^, i.e., g{k) ^ — ln|fc| in mo¬ 
mentum space. Thus we assume that the dominant con¬ 
tribution in Eq. are pairwise correlations and that 
the state can be considered as a BCS correlated state of 
bosons (where the correlations are expressed in terms of 
a permanent and not the more usual determinant). 

The assumption of only pairwise correlations is by no 
means self-evident, and there are two reasons why it is 
justihed: (a) in our system there is no real charge - pseu¬ 
dospin separation or, at least we have to allow for such 
a possibility, and thus “halves” of the (fused) quas ipar- 
ticles [in Appendix [A| the quasiparticle in Eq.( |A3[ ) is a 
product of two fused “halves”] in the pseudospin sector 
are also allowed and should be captured in the critical 
boson description, and (b) the second reason is heuris¬ 
tic: by allowing the reduction to pairwise correlations we 
can describe the ground states on the torus. The pairing 
allows the constructions of “halves”, e.g. “halving” the 
Laughlin quasihole into two non-Abelions in the Moore- 
Read stateJ^, and by keeping only pairing correlations the 
degrees of freedom that would be forgotten (suppressed) 
otherwise will be released. 

Thus we assume that the low-energy states of the in¬ 
termediate phase are in one to one correspondence to the 
states of g{r) ~ l/r^ BCS paired critical bosons. Now 
we proceed similarly to the case of permanent pairing 
g{r) ~ \jz in Ref. [THl In the effective BCS description 



















the occupation number of bosons at momentum k is 


(rik) = 2 


|g(k)P 

i-|g(k)P’ 


(15) 


where we note the minus sign with respect to the 
fermionic case. To maintain the positivity of the expres¬ 
sion with g(k) = — Alnfc, where A is a constant with 
respect to k, we must have k 7 ^ 0. This means we 
must be in the antiperiodic sector on the torus along at 
least one direction {x or y). Further, we also must have 
A ~ 1/lnfcmm because kmin ~ tt/T 1/N. Finally, 
because A —>■ 0 as —>■ 00 , i.e., no pairing in the thermo¬ 
dynamic limit, we can only have Bose condensation (at 
kmin^ • 

Let US assume (as in Ref. ESI) that we are in the (—, -b) 
sector, i.e., antiperiodic in x direction and periodic in y 
direction. We assume an even number of particles (i.e., 
that we are in the even sector), and occupy (-|-7r/L, 0) 
and (—7r/L,0) with equal probability, i.e.. 


= 'J -^/2 and = \/A^/2 6k,-krmr,^ 

(16) 

where kmin = tt/L, and bka-,o' =t,i denote second- 
quantized bosonic operators. 

Thus in the neutral sector each particle would have a 
description in the form of a two-component spinor: 


^—iQx 


(17) 


where Q = -kIL. On the other hand, in Park’s states, 
Eq. (9) in Ref. [51 we have the effective description: 


^i2Qx 

1 


(18) 


that successfully reproduces the momenta that we see in 
EDs on the torus. We may notice that the difference 
comes from the necessity to multiply the neutral part in 
Eq.(17l by exp(z(5a;) in order to “glue” the charge (pe¬ 
riodic) and neutral (antiperiodic) sector, and stay in the 
periodic sector. This “gluing” effectively brings non-zero 
momentum to the state: = NQ = Nn/L. In Park’s 

picture the momentum count results from excitons turn¬ 
ing into dipoles (not locally neutral particles), and each 
exciton acquiring the momentum k = 21:jL. Thus, at 
small length scales, charge and spin are not separated, 
and this corresponds to the requirement of gluing, i.e., 
staying in the charge-periodic sector by an extra connec¬ 
tion between charge and spin sectors. At the end, we 
have a spiral winding with q = 2Q = 27r/L, i.e., the 
spins wind exactly once over the length of the system, 
L, in X direction. Thus, due to the intertwined spin and 
charge degrees of freedom at small length scales, we ob¬ 
tain ground states at large non-zero momenta. 

But there is another important relation of the “gluing” 
procedure with what we can infer from EDs on torus. 
The description of the intermediate phase that is based 


only on the neutral sector is incomplete. As usual in 
FQH systems, the degeneracy of the ground states on 
the torus can be connected with “large gauge transfor¬ 
mations”, which are flux insertions of gauge (internal) 
fields through the two holes of the torus^. This is equiv¬ 
alent to changing the boundary conditions, but also to 
the creation of particle - antiparticle pairs. These trans¬ 
formations can create large momenta that characterize 
the degenerate ground states, or alternatively different 
sectors of the theory [see, for example, a recent discus¬ 
sion in Ref. (53]. Thus we also have to take into account 
the charge part of the theory in order to describe the 
ground state. 

Since the work of Moon et am, we know that vortex 
excitations in the exciton superfluid phase are charged 
- they are merons of unit vorticity (±1) and charge 
±(l/2)e. We may then question the CFT approach for 
constructing excitations restricted to the neutral sector 
[see Appendix [A] , or in other words whether the excita¬ 
tions obtained by Eq. (A3) comprise all possible vortex 


excitations in the intermediate critical phase. 

In answering this question we may start from the model 
function for Park’s states [Eq. (8) in Ref. [5| or Eq. D3 
here]. This construction represents two merons with op¬ 
posite charge and located at two opposite boundaries 
[compare with Eq.(135) in Ref. [S]. On the disk, this 
wave function represents one meron at the center of the 
disk and the other at the infinity. (This is somewhat 
similar to the construction of the non-Abelian quasipar¬ 
ticles in the case of Moore-Read state which also come 
in pair^^.) In this way, because the meron boundary 
constructions define the ground states, we may expect 
that merons form the quasiparticle content of the criti¬ 
cal theory. Charge and spin are not separated in these 
constructions that describe the ground states of the crit¬ 
ical phase. Thus it was necessary to include the “gluing” 
when discussing the ground state from the point of view 
that treated the charge and spin sector separately. 

More simply, on the disk with a fixed number of up 
and down particles we model a meron as 


W{zt -w)^gs{z^,Z^), cr=t,i, 


(19) 


where and w denote the ground state and the 
position of the meron respectively. Now we see the 
need for the square root factors, Hi \/[zi~~^ = 

rij \J(zj - w) Hi \/{zf - w), in the charge part (“a 
charge boost”) to relate to the expression in Eq.( [A3| ), 
which represents a neutral combination of two merons. 


B. Competing ground states on the torus and 
critical bosons 

To understand better the competing (ground) states 
present in EDs of the QHB on the torus in the interme¬ 
diate region, we will consider in more detail the odd and 
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even sectors. As in the previous subsection, we will relate 
the description of states in terms of Bose condensates at 
non-zero momenta with meron-antimeron constructions. 
For further discussion concerning the validity of the de¬ 
scription by these states, see Appendix [P] 

In the odd sector, for the periodic boundary condi¬ 
tions in both directions, i.e., in the (+,+) sector, we 
have simply Bose condensation m. K = (0,0) momen¬ 
tum. The relevant construction in the odd (—, -I-) sector, 
i.e. antiperiodic in the x direction and periodic in the y 
direction, is 


N-l 


n 




■ 1 ■ 

1 

'T 

0) 

_1 


-1 

0 

_I 


( 20 ) 


where we fixed = 1 in the ground state. 

Again, we need the gluing procedure with the boost 
K = 0). The translation by —N in the ^-direction 

followed by inversion lead to the new relevant state 
K = (A^^O). Taking also the possibility of the con¬ 
struction along the y—axis, we get a 4-fold degeneracy. 
In the even sector a similar procedure will lead to only 
the 2 -fold degeneracy for the constructions along any one 
of two directions. 

Similarly we can consider the following constructions in 
the odd (—, —) sector, i.e. antiperiodic in both directions. 


N-l r 


n 

2=1 

N-l 

n 


Qifixi+Vi) 


■ 1 ■ 

Q-if{xi+yi) 

i 

0 


2=1 




■ 1 ■ 

g-if (xi-j/i) 

2 

-1 

0 

_1 


N 


J N 


( 21 ) 

( 22 ) 


and boosts are K = (A^, A^) and K = (A^, ~^~*~^ ), 
respectively. By applying the translations by ±A^ and 
inversion, we may conclude that in the odd sector we 
have a 4-fold degeneracy with K = (A^, A^). On the 
other hand, a similar construction in the even sector leads 
to a unique state with K = (A^ 

Thus we can conclude that the prediction for the num¬ 
ber of competing ground states in the odd sector, based 
on the critical boson condensates, is in an agreement with 
particular Figs. and and Eqs. 


EDs in Section IV 


(13) and (|14|), for the description of the fc 7 ^ 0 states. 


On the other hand. Fig. shows the low lying excita¬ 
tion spectrum in the even sector with TV = 12 electrons 
at distance d = 1.51 b between the layers. We may no¬ 
tice that with respect to our expectations in the previous 
considerations (based on meron-antimeron constructions 
of ground states) there are extra low-lying states. This 
will be discussed in more detail in the next subsection. 


C. Discussion: Topological characterization 

In the previous subsection we could notice that al¬ 
though the composite boson view is convenient for the 


description of the low-lying states, it is not complete 
and sufficient for their description. In connection with 
this, we may comment that the meron construction in 
Eq. (19) seems more natural if 4'^^ is given by Eq. ^ 
(fermionic representation) than by Eq. ([^ (bosonic rep¬ 
resentation). In the fermionic case the meron, with the 
deficit of a (layer) polarization, may be interpreted as 
a Laughlin quasihole excitation of the condensate of a 
particles. 


In the absence of pairing, represented by the determi¬ 
nant in the composite fermion wave function Eq.([^, we 
could consider a bosonic system described by Halperin’s 
220 stat^ with two components, and merons as half¬ 
flux quantum, spinful quasiparticles. For that system we 
know the number of the degenerate ground states on the 
torus; we can consider all possible configurations (includ¬ 
ing meron ones) between the two edges of a cylindei!^, 
and come to the conclusion that there are 4 sectors (4 
degenerate states on a torus). They can be classified as 
members of the chiral algebra, in which we have the basic 
sector that we can associate with the identity operation, 
two sectors due to two kinds of meron - (anti-meron) exci¬ 
tations, and a sector that is associated with the one-flux- 
quantum excitation i.e. Laughlin quasihole, which can be 
viewed as a fusion of the two (f and ),) different merons. 
In this way the four sectors can be viewed as a direct 
product of two subsystems, each at filling 1/2 and having 
two sectors, due to the identity and Laughlin excitation 
in each system. Symbolically, if we denote the identity by 
1, and t meron by m, and ], meron by m, we have a direct 
product (1, m) x (1, m) = {(1,1), (m, 1), (1, m), (m, m)}, 
i.e., 4 sectors. Having the extra determinant as in Eq.(|^, 
beside the (220) factor, will not change the type of the 
charged excitations that exist in the intermediate phase, 
i.e., we expect merons. Thus we obtain a total 4-fold 
degeneracy of the ground states on the torus. 

But we have to be aware that our system is far from be¬ 
ing a simple bosonic (220) state. There we could consider 
excitations of a single boson, f or but this does not cre¬ 
ate an extra sector in the (220) system because the two 
subsystems, f and I, are uncorrelated, and everything 
boils down to the direct product of f and | (sub)system. 
In the case of the state in Eq.([^, because of the interlayer 
correlations, we may assume that a single fermion exci¬ 
tation, which we will denote by / for f subsystem and / 
for 4 , subsystem, is a non-trivial excitation, and may de¬ 
fine a sector. Thus we may consider (1, m, /) x (1, fh, f) 
as a way to generate all sectorJi^. But we notice that 
(/, /) is equivalent to two (m, ffi) excitations, and should 
not represent an extra sector. Thus, in this way we ex¬ 
pect an 8-fold degeneracy in the even sector, and this is 
corroborated by the results in Fig. 

This is by no means a usual characterization of a topo¬ 
logical phase. The topological characterization is made 
difficult by the absence of a simple model Hamiltonian 
for the model wave function [Eq. (|^ or Eq. (|^]. Due 
to (topologically) highly intertwined (pseudo)spin and 
charge, we might expect a further reduction of the num- 
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ber of relevant states. Namely, an operation of equiva¬ 
lence may be introduced for two states that are connected 
by a gapless branch, like in Fig. Nevertheless, we find 
that the low-lying spectrum, in both even and odd sector 
is robust and independent of system size, up to TV = 14 
particles in the even sector, and up to iV = 15 in the 
odd sector, and it is not in the accordance with the naive 
expectation of 4 degenerate ground states in the even 
sector that can be found in the previous literature, or as 
follows from the charge-pseudospin separated picture of 
the intermediate state. 

By comparing the tentative degeneracies on the torus 
(8 states in the even sector and 9 states in the odd sector) 
with the predictions in the even and odd sector in sec¬ 
tion VB on the basis of g{r) ~ 1/r^ BCS paired critical 
bosons, we can conclude that this view is partially suc¬ 
cessful. We may expect Abelian merons and composite 
fermions in the quasiparticle spectrum. 


also there is a drastic change in the calculated slope (ve¬ 
locity) of Goldstone mode according to Ref. [TH but the 
first excited state is smoothly connected to the Goldstone 
branch for any d. 

In contrast, on the torus [see, e.g.. Fig. [^a)] there are 
level crossings in the low-lying excited spectrum around 
d ^ is- New states appear with lower energies than the 
first excited Goldstone-branch state. These states consti¬ 
tute candidate ground states of the intermediate phase. 
They are defined by merons, as explained previously in 
Sec. |VB| Thus merons belong to the quasiparticle con¬ 
tent of the new state. 

The intermediate state has the interpretation in terms 
of BGS pairing of composite fermions (Ref. H]), and thus 
similarly to the Moore-Read case we expect half-flux 
quantum (one flux quantum in the language of super¬ 
conductivity) excitations. 


D. Discussion: Meron deconflnement 

In this subsection we summarize the reasons why 
meron deconfinement is expected in the intermediate 
state in the QHB. By “meron deconfinement” we mean 
the phenomenon that a meron can exist as a single, 
gapped excitation of the half-flux quantum and charge 
(l/2)e in the intermediate phase (as opposed to being in 
a bound pair). 



FIG. 6: The evolution of the energy spectrum in the odd 
sector on the sphere as a function of bilayer distance. Data 
is for A = 9,11,13 electrons in total (left to right), for fixed 
shift 5 = — 1 corresponding to the 111 state. 

When we compare the evolution of the QHB with 
respect to distance, on the sphere [Fig. and on 
torus [Figs. [^ [^, we see that on the sphere the low¬ 
est Goldstone-branch state evolves without level crossing 
with distance. The calculations on the sphere in Fig. 
are performed at the fixed shift 5 = — 1, which is the 
same as that of the 111 state. Note, however, that the 
decoupled Fermi liquids have a shift of 5 = —2, therefore 
the sphere geometry cannot directly capture the transi¬ 
tion without introducing some bias towards one of the 
phases. In Fig. the ground state smoothly changes 
around d £b due to the change in the correlations, and 


VI. CONCLUSIONS AND DISCUSSION 

We presented a description of the intermediate region 
of the QHB in terms of critical bosons based on detailed 
exact diagonalizations on the torus. Our study comple¬ 
ments earlier results on the spher^ and torufl geometry, 
where an intemediate phase was identified at bilayer dis¬ 
tances d « 1.5£b- These studies have pointed out that 
the phase involves pairing of composite fermions, and has 
pseudospin spiral order. Our work shows that this state 
can also be related to critical bosons and is likely gap¬ 
less. The critical condensates at finite wave vectors of 
the intermediate phase on the torus can be related to the 
existence of merons as the lowest lying charged quasi¬ 
particles, i.e., to meron deconfinement. The topologi¬ 
cal ordering in the intermediate region with meron frac- 
tionalization allows a smooth transition (ordering) in the 
adjusted phases in which merons are confined in pairs; 
as vortex excitations in the XY ferromagnet for small 
d, and, for large d, as bound pairs they make neutral 
fermions. 

We caution that the critical boson representation em¬ 
ployed here may be less powerful than similar descrip¬ 
tions of the more conventional FQH states where charge 
and pseudospin (neutral degrees) are not so intertwined. 
In our case, the natural representations in terms of 
(composite) bosons (exciton superfluid) and (composite) 
fermions (composite Fermi liquids) hold in the limiting 
cases of small and large distances. However, we showed 
that most of the low-energy features of the intermedi¬ 
ate state can be conveniently explained in the picture of 
critical bosons. 

Assuming the intermediate state develops into a crit¬ 
ical point in the thermodynamic limit, the QHB system 
may offer a realistic application of the ideas of deconfined 
criticalitj^^. Moreover, the presented QHB description 
may also explain the phenomenon of “meron loosening” 
in experiments^. Before invoking any disorder effects, 
we may assume that merons are liberated (not bound 
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in pairs with opposite layer indexes), and equally con¬ 
tribute to inter and intra dissipative processes that lead 
to the breakdown of the quantum Hall effect. As argued 
in Ref. HTJ the quasiparticles with one-half of the elemen¬ 
tary charge are natural (confined in dipoles) constituents 
of the underlying composite fermions in the Fermi-liquid¬ 
like state (at filling one-half), due to the particle-hole 
symmetry. In our case these are merons (of the same 
layer index), which are brought to the existence in the 
intermediate state, before the two separate Fermi-liquid¬ 
like states (each at filling one-half) set in. 


Appendix A: The bosonic interpretation 


We can 




rewrite the neutral part in Eq. 


0 


as 



Thus, the neutral part in Eq. 0 in the CFT language is 
the correlator of equal numbers of two vertex operators: 
gZ0(^)-Hi0(2) g-*0(2)-i0(2) denote 

the holomorphic and anti-holomorphic part of the (non- 
chiral) bosonic field, cj){z,z) = (j){z) + (j){z). The basic 
correlator is: 


^gi/3(^(zi,2i)g-i/30(z2.Z2)^ 


1 


(Al) 


We expect that the excitations are also represented by 
vertex operators. For their construction we have two pos¬ 
sibilities: and and the corre¬ 

sponding conjugates, which in charge-neutral combina¬ 
tions in the correlators with basic particles 
and i.e. up and down (neutral) excitons, 

contribute factors. 


^i6i(f>{w,w) 




(A2) 


and, with the field, 9{z, z) = (j){z) — i4>{z), we have 


.S^eiw^w) . nP4 - ^ (.nx 

To ensure the single-valuedness of the ensuing electronic 
wave function we must have S 2 = 1/2. The same ar¬ 
gument would allow any rational including (5i = 0. 
Therefore the excitation spectrum contains a branch con¬ 
tinuously connected with the ground state and we expect 
that the state in 0 supports a branch of gapless exci¬ 
tations. (If we allow only (5i = 1/2 we would have the 
quasiparticle content of the deconfined phase of Z2 gauge 
theory.) Thus the state (|^ is a critical state with the 
quasi-long-range order. 


hm m-{zi,z,)'^7{z2,Z2)) = 


|zi-Z2|->'00 


^gi0(zi,Zl)g-i0(z2,Z2)^ _ 


Zl - Z2 


(A4) 


among up and down excitons. 


Appendix B: Projection to the LLL 


In this Appendix we explain the projection to the 
LLiP that defines the LLL projected wave function of 
p-wave paired composite fermions in Eq. 0. First we 
expand and rewrite Eq.Q as 


4- 


PSF — 


VlllI E 


(tGS 


N/2 


4'220(4)--->2^Ar/2>^' 


(^1 ^o'(l)) ■ ■ ■ (^Af/2 ^(T(Ar/2)) 

(Bl) 


^N/2 


)}■ 


where by 4*220 we denoted the part with the Laughlin- 
Jastrow correlations. 


N/2 N/2 

4-220 =- 4)" n(4 - 4)"- (b2) 

i<j i<j 


To understand the projection, we consider an arbitrary 
permutation a under the sum in Eq. (Bl). For the fixed 
permutation, we can express the product 4*220 as a poly¬ 
nomial in {{zj - {zj + ^^(i)); I = 1,..., N/2} vari¬ 

ables. In this way, when we combine 4*220 with the pair¬ 
ing part with fixed cr, the projection reduces to the pro¬ 
jection of each ^ (with m = 0,..., N — 2), to the LLL. 
If we choose the symmetric gauge, the basis of the LLL 

k 

is given by 4*fc(z) = 7=^, with k = 0,.. .,N-l. Here 

we dropped the Gaussian factor, exp{—|zp/4}, which is 
included in the measure defining the scalar productP^. 
Thus, 


,z^, , 27rTO!2"- 

Z ' V27rfc!2'= 


(B3) 


It follows that 




^m+1 


2{m + 1)' 


(B4) 


In particular, for m = 0, Vlll{\) = f, he., the leading 
term in the short distance approximation reproduces the 
behavior of the p-wave paired composite fermion wave 
function described by Eq.([^. 


Appendix C: Relation to pseudospin spiral order 

The first exact diagonalization study on the torus that 
pointed out to the existence of the competing ground 
states with diverging momenta in the thermodynamic 
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limit was done by Park in Ref. dl Our numerical results 
are in agreement with Ref. [21 in particular for iV = 13 
electrons and bilayer distance d = I.BIb [Fig. [^a)]. 

Ref. 121 pointed out that one energy level (the lowest 
excitation at non-zero momentum, which represents four 
states due to inversion symmetry) is “ripped out” from 
the quasi-continuum of high momentum states. The mo¬ 
menta that correspond to these states are the ones in Eq. 


(131 in the main text, or 


(iV^,0), {N^,0), (0,iVt)and(0,iV^), 


(Cl) 


in units of 27r/L where L = y^2Tr£'^gN. 

In order to motivate and explain the emergence of the 
particular momenta above, Park discussed the “spiraling 
states” in the Landau gauge, 

= n + cLV)|0)- (C2) 

-n » 


Here denotes the electron creation operator, with in¬ 
teger m = 0,..., iV — 1 labeling the linear momentum 
along the cylinder, i.e., 2^. With respect to the ground 


state of the XY ferromagnet [n = 0 in Eq. (C2), or Eq. 
Q if we consider the fixed number of particles represen¬ 


tation], the acquired momentum is = ( 


27rnNt 


, 0 ). 


In the first quantization, in terms of the usual Landau 
gauge states 


exp 


/ 2'!rimzk \ 

I L j 


exp 


1 


(C3) 


where Zf^ = Xk + k = 1,..., N, we can express the 
spiraling state as 


n 


27r. 


exp (j ^ 
1 


P[(Zi-Z^)exp{-i^y2}, (C4) 


i<j 


where we introduced Zj = exp ^ ^ Thus the physics 

of the neutral sector is essentially given by the factor 


n 


exp 

1 


(C5) 


which describes long spiraling with period L in the pseu¬ 
dospin space. 

Park argued that by more elaborate inclusions of the 
disordering effects in the ground state, similar construc¬ 
tions to the ones in (C2) should be relevant for the de¬ 


scription of the new ground state(s). They will corre¬ 
spond to the new intermediate phase with the “pseu¬ 
dospin spiral order”. 

In this work we point out that there are extra states, 
as seen in Pig. falling down, which also belong to the 
set of candidate ground states of the intermediate phase. 
Given the unusually large momenta (and long periods) 
for spiral order, we argue in Sec. |V] that they are very 
likely an expression of the topological ordering in the in¬ 
termediate phase. Furthermore, we connect them with 


merons - the spin-texture fractional excitations of the 
XY ordered phase - which are confined (bound in pairs). 
The long spiraling in the ground states of the intermedi¬ 
ate phase on the torus is thus naturally explained by the 
spin-texture description of merons. 


Appendix D: Description of the low-lying modes 

To find and understand the description of the low-lying 
modes it is convenient to begin with the case d = 0 
{SU{2) symmetric case) in the presence of only Vq pseudo 
potential for which we can identify the whole family of 
zero modes (zero energy states). As shown in Ref. [21] the 
basic states (from which we can get all other zero-energy 
states) similarly to edge states in the Laughlin case, can 
be described in the disc geometry as 



T t 

Z- Z- 




(Dl) 


i<j 


where the sum is over various L—tuples of f particles. 

With the inclusion of higher angular momentum 
pseudo potentials they will not be degenerate anymore, 
but will constitute the candidate states for the states 
of the lowest energy at a given momentum L (with re¬ 
spect to the ground state momentum). For example, the 
state at L = 1, ~ -^i)) can be associ¬ 

ated with the first excited state above ground state in 
the Goldstone branch. This was demonstrated in Ref. im 
when a boost of angular momentum L acting only on one 
pseudo-spin projection (layer index) in the ground state 
was identified with the Goldstone mode with momentum 
L. This ansatz was very satisfactory for the whole range 
of distances (from zero to infinity) in the bilayer for low 
momenta (L = I, 2, 3). 

On the other hand the state in Eq. (Dl I when L = 


Nf 


i<j 


(D2) 


i=l 


is a meron-(anti-meron at the infinity) state. This state 
is a projection of the meron state introduced in Ref. [5] 




I „Tt 

n ' ^n+1 


)|o), 


(D3) 


where =t, j. are electron creation operators in the 

second-quantized notation, to the fixed relative num¬ 
ber, of electrons. This state was also used by 

ParlP to model the state(s) that fall(s) from the higher- 
momentum region when the distance (d) increases. 

The best of overlap in Ref. |2] was reached around d < 
£b - after that the overlap was rapidly decreasing. As the 
author already pointed out, the reason for this is that the 
used construction is based on the 111 state, i.e. the state 
in Eq. [^projected to the fixed number of electrons, and 
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we need to include the changes in the ground state that 
set in around d ^ £b- Thus we may imitate the approach 
of Ref. [in] and model new state and correlations by 


(D4) 

i=l 


where 'I'o represents the exact ground state at distance 
d that we find in fc = 0 sector. 

Nevertheless, the quantum numbers of relevant states, 
Eq. D4 are the same as from the states they evolve from, 
Eq. D2 and we can use these states i.e. their description 
given by Eq. [^to classify all possible low lying states on 
the torus. 


Appendix E: Critical region 


Our analysis suggests that the bilayer intermediate 
state appears as a critical region between a ferromagnet- 
ically ordered state and a (completely) disordered state. 
It is tempting to expect that by modifying the param¬ 
eters in the critical region we can reach the Z 2 D phase 
(the deconfined phase of Z2 gauge theory, mentioned in 
Sec.|VD|. This is one of the simplest topological phases 
that would be realized by strong pairing of composite 
bosons. Also, it is important to understand better the 
influence of the various Haldane pseudopotentials on the 
competing ground states in the critical region. 

In order to explore the phase diagram more broadly, 
we consider modifications of the bare Coulomb potential 
by the short-range Haldane pseudopotential term^ 




2{2n)Y,SV^Lm{q^), (El) 

m 


where m = 0,1, 2,.... The purpose of these modifications 
is to model the perturbations that occur in an actual 


experimental system (for example, effects due to finite 
width of the sample in the perpendicular direction, mix¬ 
ing with higher Landau levels, etc.). Note that we are 
free to separately tune the intra- and inter-layer interac¬ 
tions. Adding even m pseudopotentials will only have an 
effect on inter-layer potential because the Fermi statistics 
inside each layer precludes even values of m. 

In Fig. we can see the evolution of the spectrum 
with small departures {SVm\rn = 0,...,5} of the Hal¬ 
dane pseudopotentials from their Coulomb values. The 
influence of the even pseudopotentials, (5Vb, i5y2, and SV 4 
is qualitatively the same; for stronger positive values the 
ferromagnetically ordered state is stabilized, and for neg¬ 
ative values we expect a compressible behavior. The 
influence of 5Vi is according to our expectations; for 
stronger positive 6 Vi we see compressible behavior, while 
for negative 5Vi the ordered state is stabilized. 

On the other hand the influence of positive 6 V 3 and 
6 V 5 is interesting as it stabilizes one of the set of states 
(competing ground states) that have characteristic mo¬ 
menta that are compatible with Z 2 D. Nevertheless, a 
closer inspection of the influence of 5Vz does not lead to 
a firm (topological) characterization of the ground states 
for positive values of 5V^. Moreover, the behavior of the 
gap in this region is strongly dependent on system size, 
which further suggests the phase is compressible. 
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FIG. 7: (Color Online) Energy spectrum for the Coulomb interaction at bilayer distance d = 1.5£b modified by short-range 
Haldane pseudopotentials 5Vm, m — 0,... 5 [(a)-(f)]. In all cases, total number of electrons is = 11 and the unit cell is 
a square. The spectrum is obtained by finding 10 lowest eigenvalues in each momentum sector, and for clarity purposes the 
resulting energies are plotted relative to the ground state for each value of the pseudopotential. Levels belonging to high- 
symmetry momentum sectors are denoted by special colors. For odd pseudopotentials 5Vi,5V2.,5V^, we modify the intra-layer 
and inter-layer Coulomb by equal amounts. The even pseudopotentials 5Vo, 514, 514 only affect the inter-layer Coulomb because 
of Fermi statistics. 
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